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H(x )  heaviside  step  function 


JtUA)  kernel  function  (see  Eq.  (2.2)) 

conventional  slender  body  kernel  function,  see  Eq.  (2.24) 
K(k)  elliptic  integral  of  first  kind  (see  Eq.  (2.5)) 


see  Eq.  (2.6) 


total  lift,  see  Eq.  (2.10a) 
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L(x) 


normalized  load  function 
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R'(x) 
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Tn(x) 

Un(x) 
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x,r,0 


load  function  (see  Eq.  2.9)) 

^  -  Xjj,  body  length 
normal  force  per  unit  length 

total  moment  about  center  of  body,  see  Eq.  (2.10b) 

| vQ|  surface  speed 

reference  body  radius 

body  radius,  see  Fig.  2.1 

dR/dx  ,  body  slope 
2 

,  body  cross -sectional  area 
Chebyshev  polynomial  of  the  first  kind 
Chebyshev  polynomial  of  the  second  kind 
body  volume 

axial  flow  surface  velocity 

free  stream  velocity 

upwash  function  (see  Eq.  (2.8)) 

axial  locations  of  the  nose  and  tail  of  the  body  of  revolution 
see  Eqs.  (3.28)  and  (3.29) 
cylindrical  coordinate  (see  Fig.  2.1) 


a  angle  of  attack 

c  2Re/t 

p„  free  stream  density 
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2R(x)/t,  normalised  body  radius 


I.  INTRODUCTION 


The  viscous  theory  of  lift  on  bodies  of  revolution  is  developed  in  detail 
in  Ref.  1.  The  problem  of  calculating  lift  on  a  body  is  reduced  to  one  of 
solving  an  integral  equation  of  the  first  kind,  in  complete  analogy  with  the 
theory  of  lift  on  wing  sections.  Whereas  the  theory  is  formally  well 
developed,  specific  results  are  difficult  to  obtain  due  to  the  complexity  of 
the  associated  kernel  function.  In  Ref.  1,  an  explicit  evaluation  of  the 
kernel  function  was  carried  out  and  load  calculations  were  made  for  the  case 
of  a  ring^wing,  a  degenerate  but  highly  useful  form  of  a  body  of  revolution. 
The  results  obtained  were  in  agreement  with  a  theoretical  result  of  Ribner 
(Ref.  2)  and  with  experimental  data  of  Flatau  (Ref.  3). 

In  the  present  report,  the  theory  is  simplified  by  deriving  the 
asymptotic  form  of  the  kernel  function  for  slender  bodies  of  revolution.  The 
resulting  inviscid  integral  equation  is  solved  in  closed  form  for  the  case  of 
a  prolate  spheroidal  body.  An  explicit  representation  of  the  inviscid 
eigenfunction  is  obtained  and  numerical  results  for  the  lift  distribution  are 
compared  with  results  of  conventional  slender  body  theory.  Also,  the 
magnitude  of  the  eigenfunction  is  estimated  by  correlating  with  experimental 
data  on  the  USS  Akron.  The  qualitative  features  of  the  load  distribution  are 
in  much  better  agreement  with  experimental  data  than  the  distribution  obtained 
with  the  equivalent  base  area  concept  applied  to  conventional  slender  body 
theory. 


t 


II.  ASYMPTOTIC  STEADY  STATE  THEORY 

The  point  of  departure  for  the  present  study  is  the  inviscid  integral  Eg. 
(2.35)  of  Ref.  1;  i.e.. 
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^g(x)/V. 


(2.1  ) 


where  )  is  the  unknown  load  function, 

function  and  )  is  the  kernel  function.  The 

state  kernel  function  is 


is  a  specified  upwash 
definition  of  the  steady 


where 


j^0c,5  ) 


Lim 
r+R(x ) 


R(5'  ) 
P=R(£') 


ae  • 


cos0  d0 


/x2  +  r2  +  p2  -  2rp  cos0 


(2.2) 


JTl 


x  +  (r*p )‘ 


(-*) 


K(k)  +  —  E(k) 


with  the  standard  elliptic  integrals  (see  Ref.  4,  p.  589) 


E(k) 


/*/2  - 

/l  -  k2  sin2© 


d6 


(2.3) 


(2.4) 


II 


(2.5) 


where 


r 


k2  _  **~P 

x2  +  (r-K>)2 


(2.6) 


Also 


3  3  3  3 

—  -  -  R'(x)  —  +  ,  ~ 

3N  dx  3r  3N> 


3  3 

-“'“''Ir  *57 


(2.7) 


and  R(x)  is  the  specified  body  and  wake  control  surface  shape  (see  Figure 

2.1). 


For  the  remainder  of  this  report  we  consider  a  body  of  revolution  at 
constant  angle  of  attack  a.  To  lowest  order  in  body  thickness  the  upwash 
function  becomes 


(q0/v.) 
/l  +  R'2 


(2.8a) 


«  -  a  for  slender  bodies 


(2.8b) 


The  approximation  made  in  Eg.  (2.8b)  is  not  valid  near  the  fore  and  aft  ends 
of  the  body  where  the  flow  stagnates.  In  these  regions  one  should  use 
Eg.  (2.8a)  with  the  calculated  or  measured  axial  flow  surface  speed,  gQ. 

The  unknown  load  function  J^*fx)  in  Eg.  (2.1)  is  related  to  the  normal 
force  per  unit  length  on  the  body  l(x)  via  the  following  relation: 


.  ,  1  2 
*(*)  »  2  P“V"  * 


( 2nR(x) )  .^(x) 


(2.9) 


The  total  life  and  moment  about  the  y  axis  are  given  respectively  by 
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l(x)dx 


(2.10a) 


and 
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Control 

Surface 


r 
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x* A (x)dx 


(2.10b) 


Let  A  *=  (xT  -  xN)  denote  the  total  body  length  and  define  the  normal  force  and 
moment  coefficients  as  follows: 


_ L _ 

1/2  P«,vJ*xR0*A 


Normal  Force  Coefficient 


CM  =  - - -  Moment  Coefficient  (2.11) 

1/2  P»v21tRo.a.(A/2) 

where  Rq  is  a  reference  body  radius;  e.g.  the  radius  at  the  maximum  girth  or 
center  of  buoyancy.  Substitute  Eq.  (2.9)  into  the  formulas  (2.11)  to  obtain 

and 

CM  =  *  j  Jz^x>  *  x  *  *  dx  (2.12) 

In  the  last  two  formulas  and  in  the  remainder  of  this  report  we  choose  the 
body  half  length  (A/2)  as  a  reference  length  and  further  choose  the  origin  of 
coordinates  at  the  body  mid-length.  All  integrals  over  the  body  then  become 
integrals  from  -1  to  1  in  normalized  coordinates.  The  moment  is  about  the 
center  of  the  body  and  is  positive  nose  up. 

The  most  important  step  in  reducing  the  integral  Bq.  (2.1)  to  a  more 
tractable  form  is  to  introduce  the  slender  body  approximation;  i.e.r 

R(x)  <<  A/2  (2.13) 
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The  kernel  function  (see  Bq.  (2.2))  can  then  be  approximated  by  a  relatively 
simple  expression.  For  load  points,  C,  that  are  far  removed  from  the  upwash 
control  point,  x;  i.e.. 
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|x-d  >>  R(x ) 


(2.t4) 


the  kernel  function  can  be  replaced  by  the  classical  slender  body  kernel  as  we 
demonstrate  below.  For  load  points  close  to  the  upwash  control  point;  i.e., 

|x-e|  «  R(x)  (2.15) 


the  kernel  can  be  replaced  by  a  local  ring-wing  approximation.  We  turn  now  to 
a  detailed  derivation  of  these  results. 

For  the  case  covered  by  Eq.  (2.14),  we  recognize  that  the  parameter  k  in 
Eq.  (2.3)  is  much  less  than  unity  so  that 


^(k.r.p) 


k  «  1 


(2.16) 


and 


R(5') 


32^  irR(g’) 
3N3N*  "  | x_£ *|  3 


(2.17) 
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Thus,  for 


x  «  5 


we  have 


«^(x,£ ) 


d£  ' 

(t'-x)3 


X 


a  (2.18) 

2(£-x)2 


ft 
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For 


x  >>  £ 


vie  write  Eq.  (2.2)  in  the  following  form: 


Lim 

r+R(x) 


3N3N ' 


R(C')  d£* 


■fJ0L^ 


*  f(x) 


ttR(£  ) 

2 (x-{ ) 2 


(2.19) 


where  the  second  term  is  obtained  with  the  same  approximation  that  led  to 
Eq.  (2.18).  The  function  f(x)  can  be  inferred  from  our  ring-wing  results  as 
follows.  Note  that 


f(x) 


Lim 

r+R(x) 


jm: 

3n3n  ' 


R(C’)df 


(2.20) 


and  furthermore,  that  the  principal  contribution  to  the  integrand  is  for 
points  near  x.  Thus,  we  replace  R(S')  by  R(x)  and  let 


_af _ §1 

3N3N '  “  3r3p 


in  Eg.  (2.20).  The  result  is 


f(x)  -  Urn  R(x )  f~  d£' 

r,p+R(x)  3r3p 


(2.22) 


_ *_ 

R(x) 


(2.23) 
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where  we  have  noted  that  Eq.  (2*22)  is  the  asymptotic  form  of  the  ring -wing 
kernel  (see  Eqs.  (2.42)  and  (2.45)  of  Ref.  1). 


With  Eqs.  (2.18)  and  (2.19)  we  can  express  the  kernel  function  of 
classical  slender  body  theory  in  the  form 


(2.24) 


If  we  neglect  the  second  term,  the  integral  equation  of  slender  body  theory 
becomes  (here  we  use  dimensional  coordinates) 


3&H 


or 


so  that 


Ruxyv,) 

A  +R'2 


(2. 25) 


^^(x) 


d  (R(x)(qo/v,)  ) 
dX  (  A  ♦  R’2  ( 


(2.26) 


Substitute  the  last  result  into  Bq.  (2.9)  to  obtain  the  lift  per  unit  length. 


*o(x)  "  p«vi  *  * 


R(x)(qo/v.) 


d  lRlit>(yvJ 

dx  A  +R'2 


(2.27) 


If  all  second  order  effects  of  slope  are  neglected,  then  we  obtain  the 
classical  slender  body  formula. 


tQ(x) 
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dS 

dx 


a 


(2.28) 
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with 


S  *  ftR2(x)  (2.29) 

A  slight  modification  of  Eq.  (2.28)  that  forces  the  lift  per  unit  length  to 
zero  at  blunt  ends  is  obtained  by  replacing  the  body  radius  R(x)  by 
R(x)(^0/v«).  This  is  tantamount  to  neglecting  the  factor  / 1  +  R'2  in  the 
denominator  of  Eq.  (2.27).  Alternatively,  we  can  replace  the  dynamic  pressure 
in  Eq.  (2.28)  by  its  local  value  Since  we  do  not  rigorously 

calculate  second  order  corrections  for  body  slope,  all  of  the  foregoing 
approximations  are  formally  equivalent. 

We  pointed  out  in  Ref.  1  (p.  17)  that  the  most  significant  nonuniformity 
of  classical  slender  body  theory  is  due  to  the  ad  hoc  elimination  of  the 
Cauchy  singularity  in  the  fundamental  integral  Eq.  (2.1).  (See  Figure  2.2  for 
a  comparison  of  singularities  in  the  various  kernels. )  The  main  theoretical 
result  of  this  study  is  the  correction  of  the  classical  slender  body  kernel 
(Eq.  (2.24))  that  retains  the  essential  Cauchy  singularity.  This  step 
permits  us  to  calculate  the  inviscid  eigenfunction  and  to  introduce  the  effect 
of  viscosity  to  determine  the  magnitude  of  the  eigenfunction. 

To  obtain  the  correct  lowest  order  slender  body  kernel  function  we 
evaluate  Eq.  (2.2)  for  | x— Cl  <<  R(x).  We  have  already  seen  from  our  ring-wing 
calculation  that,^fx,£)  is  singular.  Furthermore,  we  can  use  the  ring-wing 
result  directly  since  to  lowest  order  in  thickness,  a  slender  body  can  be 
approximated  locally  be  a  cylindrical  surface.  Thus,  from  Eqs.  (2.41)  and 
(2.46)  of  Ref.  1, 


^(x,C)  * - -  for  lx-€|  «  R(x) 

x-i 


(2. 30) 


The  complete  asymptotic  kernel  function  is  obtained  by  adding  Eq.  (2.30)  and 
the  principal  part  of  Eq.  (2.24).  Thus 


2n 


B"-?1 


(2.31 ) 


Substitute  Eq.  (2.31)  into  Eq.  (2.1)  and  multiply  the  resulting  equation  by 
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(-R(x))  to  get  (in  normalized  coordinates) 


T  (x) 
n 


L(y) 

x-y 


dy  + 


L(y )dy  =  x (x)  a (x) 


(2.32) 


where 


Vv« 

O(x)  *=  "  .  ,  :=  •  O  “  O 

/l  +  R'2 


(2.33) 


T  (X  ) 


2R(x) 

£ 


(2.34) 


and  L(x)  is  the  load  function  Jz(x )  expressed  in  terms  of  the  dimensionless 
axial  coordinate. 

The  slender  body  integral  Equation  (2.32)  is  the  fundamental  theoretical 
result  of  this  report.  We  investigate  the  solution  of  this  equation  in  the 
following  section.  In  conclusion,  we  remark  that  for  x(x)  <<  1,  and  with  the 
classical  slender  body  approximation,  we  recover  the  result  of  Equation 
(2.25).  On  the  other  hand,  for  T(x)  >>  1,  we  obtain 


1  f  MU  dy  „  a(x) 

"  7-1  x-y 


(2.35) 


the  integral  equation  of  wing  theory  or  actually  the  limiting  case  of 
ring-wing  theory  with  an  infinite  diameter  to  length  ratio.  It  is  essential 
that  both  terms  be  retained  in  Eq.  (2.32)  to  obtain  correct  slender  body 
results. 
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III.  SOLUTION  OF  THE  FUNDAMENTAL  INTEGRAL  EQUATION 
FOR  A  PROLATE  SPHEROIDAL  BODY 


Consider  the  body  shape 

t(x)  *  e/l  -  x2  ,  |x|  <  1  ,  e  «  1  (3.1) 

and  suppose  that  the  wake  control  surface  (see  Figure  2.1)  has  been  collapsed 
onto  the  line  extending  downstream  to  infinity  from  the  aft  most  position  on 
the  body.  Next,  we  assume  a  load  function  of  the  form 


L(x) 


Vn<*> 

n=0 


(3.2) 


where 


T  (cos0)  *  cos  n  0  ,  x  *=  cos0  (3.3) 

n 

The  last  expression  defines  the  Chebyshev  polynomials  of  the  first  kind,  Tr(x) 
(Ref.  4,  p.  776).  Now  substitute  Eqs.  (3.1)  and  (3.2)  into  Eq.  (2.32).  The 
result  is 


e  a 


(3.4) 


where 


Un(cos0  ) 


sin(n+1 )9 
Bind 


(3.5) 


The  last  expression  defines  the  Chebyshev  polynomials  of  the  second  kind, 
Un(x).  Finally,  multiple  Eq.  (3.4)  by  UB(x)  and  integrate  over  x  to  obtain 
the  following i 

2AQ  -  (1  +  e  )A,  »  e  a  n-1  (3.6a) 
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2Aq  -  (i  +  en)^  *  0 


n-2, 3, .  •  • 


(3.6b) 


Thus,  we  obtain  for  the  coefficients  AR  that  define  the  load  function. 


2  ea 

A<  a  -  A_  —  — 

1  1  +  e  °  1  +  e 


(3.6) 


A  =  - - -  A 

n  1  +  en  ° 


n-2 1 3 , .  •  . 


(3.7) 


The  coefficient  AQ  is  completely  arbitrary. 

With  the  foregoing  results  we  can  express  the  load  function  as  follows: 


_£2_  .  .  *  +  ,F(x)- 
1  +e  ° /rr^ 


(3.8) 


where 


F(x)  -  2f(x)  -  1 


(3.9) 


£  TT 


T_  (x ) 


(3.10) 


The  first  term  in  Eq.  (3.8)  is  the  particular  solution  of  Eq.  (2.32)  and  is 
essentially  the  same  result  that  one  would  obtain  with  classical  slender  body 
theory.  The  second  term  is  the  important  new  result.  It  is  the  eigenfunction 
multiplied  by  the  arbitrary  coefficient  AQ.  We  now  derive  an  explicit 
functional  relation  for  the  eigenfunction  F(x). 

Use  the  definition  (3.3)  in  Bq.  (3.11)  to  obtain 


g(9  )  “  f  (cos8 )  ■  Re 


j?r.« 


(3.11) 
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where 
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z) 
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1  +  ne 


(3.12) 


Now  it  can  readily  be  shown  that 
equation 


satisfies  the  first  order  differential 
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zz^'(z)  +  -  — - — 

1  -  z 

with  the  boundary  condition 

Jlo)  -  1 

Thus,  we  integrate  Eq.  (3.14)  to  obtain 


(3.13) 


(3.14) 


t 
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t 


t1/e-1 
1  -  zt 


dt 


(3.15) 


From  Ref.  4,  p.  556,  we  find  that  F(z)  can  be  expressed  in  terms  of  the 
standard  hypergeometric  function;  i.e.. 


Jlz)  -  F(  1 ,  1/e;  1  +  1/e,  z) 


(3.16) 


or  the  incomplete  Beta  function  (see  Ref.  4,  p.  263) 


— —  Bt(1/e,  o) 
£ 


(3.17) 


With  either  Bq.  (3.16)  or  (3.17)  the  eigenfunction  is  completely  defined  for 
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arbitrary  fineness  ratio  (1/e).  For  the  special  case  that  e  is  the  reciprocal 
of  an  integer;  i.e.. 


e  ■  1/N  (3.18) 

where  N  is  the  body  fineness  ratio/  we  can  express  f(x)  in  terns  of  elementary 
functions.  The  result  is 


f  (cos©  )  *  -  N 


(6  -ir ) 
2 


sinN0  +  cosNS 


l n(2  sin0/2)  + 


cosne 

(N-n) 


(3.19) 


Note  that  f(cos0)  has  a  logarithmic  singularity  at  0  *  0;  i.e.,  at  the  body 
aft  end.  Also,  we  remark  that  even  though  the  eigenfunction  is  expressed  as  a 
finite  sum  of  oscillatory  trigonometric  functions,  the  actual  eigenfunction  is 
a  smooth  function  of  x  as  we  show  below. 

To  complete  the  inviscid  theoretical  work  in  this  report,  we  derive 
expressions  for  the  load  distribution  and  the  total  normal  force  and  moment 
coefficients.  Define  the  normal  force  coefficient  per  unit  normalized  length 
as 


cN(x) 


L(x) 


.  R(x) 
v-  Ro 


-  L(x) 


#  T (x ) 
v-  To 


(3. 20) 


where 


T 


O 


(3.21) 


Then  from  Bqs.  (2.12),  we  have 


§ 


t 


°N  * 


CM  *  - 


Cjj(x)  dx 


cN(x)  •  x  dx 


(3.22) 


Substitute  Bqs.  (3.1)  and  (3.8)  into  Bq.  (3.20)  and  let  tq  *  e  to  obtain 

'h'"  -  -  irfn  ■  (%r)  ♦*<>•(;:  ,,xj  ,3-23> 

Since  qQ(x)  is  a  symmetric  function  of  x  for  a  symmetrical  body,  we  further 
obtain 


CN  “  Ac 


£i 


F(x )  dx 


and 


f  f  i 

'  ' c  J-,  v-  J. i  v- 


x  •  F(x)  dx 


(3.24) 


An  important  conclusion  can  be  made  at  once  from  the  foregoing  results. 
Hie  normal  force  coefficient  is  proportional  to  the  magnitude  of  the 
eigenfunction  in  complete  analogy  with  the  theoretical  lift  coefficient  on  a 
wing  section.  Furthermore,  this  conclusion  is  not  limited  to  a  prolate 
spheroidal  body.  Die  general  form  of  cN(x)  is  the  slender  body  distribution 
of  load  plus  the  contribution  of  the  eigenfunction.  For  a  closed  body  the 
first  part  always  integrates  to  sero.  With  the  present  formulation  it  is  not 
necessary  to  introduce  an  artificial  base  area*  Experimental  normal  force 
data  can  be  correlated  in  terns  of  AQ  and  the  complete  normal  force 
distribution  can  be  calculated  with  Bq.  (3.23). 


# 
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Examples  of  the  theoretical  normal  force  distribution  are  presented  in 
Figure  3*1.  Me  actually  plot  the  approximate  load  function 


-  x  +  A£  F(x) 


(3. 25) 


where 


(3.26) 


As  a  brief  aside,  the  surface  speed  gQ  on  a  prolate  spheroidal  body  is  known 
in  closed  form  (see  Ref.  5,  p.  141,  Eq.  4  for  the  potential  function)  and  is 
proportional  to  ( 1 /J\  +  R '  ^ ) .  We  omit  all  blunt  end  mean  flow  effects  in  the 
following  comparison  in  order  to  stress  the  differences  between  conventional 
slender  body  theory  and  the  present  theory  that  includes  the  eigenfunction. 
The  end  effects  must  definitely  be  included  as  the  theory  is  refined  into  a 
working  tool. 

For  A^  «  0,  we  obtain  the  conventional  slender  body  result  (so  labeled  in 
Figure  3.1).  The  eigenfunction  F(x)  is  plotted  for  body  fineness  ratios,  N  « 
7  and  11.  Recall  that  the  eigenfunction  can  be  evaluated  with  Eg.  (3.19)  when 
c  is  the  reciprocal  of  an  integer.  The  eigenfunction  for  N  •  7  is  added  to 
the  slender  body  load  to  obtain  the  curves  for  A£  -  1/2  and  1.  It  is  an 
interesting  and  physically  appealing  feature  of  the  theory  that  the  slender 
body  load  distribution  is  very  good  over  the  forward  portion  of  the  body.  The 
nature  of  the  eigenfunction  is  such  that  the  negative  slender  body  load  on  the 
aft  section  of  the  body  is  reduced.  In  Figure  3.2  (taken  from  Ref.  6)  the 
theoretical  and  measured  load  distribution  on  the  USS  Akron  are  compared. 
Even  though  the  results  presented  are  well  into  the  nonlinear  angle  of  attack 
region  (a  ■  12  deg),  it  is  noted  that  slender  body  or  Hunk  theory  is  not  too 
bad  over  the  forward  65%  of  the  hull.  The  deviation  between  slender  body 
theory  and  experiment  also  has  the  same  qualitative  features  as  the 
differences  indicated  in  our  theoretical  Figure  3.1. 


■ 1 


g(x)=-X  +A0F(X) 


Logarithmic 

Singularity 


Slender  Body  Theory  + 
Eigenfunction  for  N  *  7  (  €  =  '/>) 


Eigenfunction 


N  =  7 


Prolate  Ellipsoid  of  Revolution  at 


fMAL  FORCE  LOADING  n 
DYNAMIC  PRESSURE  ’  Q 


t 
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We  also  present  a  comparison  of  the  equivalent  area  concept  with  the 
present  theory  in  Figure  3.1.  If  we  calculate  the  integral  of  g(x),  i.e.,  the 
total  load,  we  get 


ii-f  9(x)dx  -  J 


F(x)dx 


(3.27) 


On  the  other  hand,  the  integral  of  the  slender  body  load  up  to  station  x  is 


g»'- 


■-/; 


1  -  X. 


xdx  = 


(3.28) 


N  N 1 

If  the  equivalent  area  is  selected  so  that  gQ  -  gQ  ,  we  obtain 


z*1 

xo  = 

;  •  aij 

1  F(x )dx 

-1 

1/2 


(3.29) 


where  integrals  of  the  eigenfunction  are  given  in  Table  3.1. 

TABLE  3.1 


N 

r1 

J  F (x)  dx 

r1 

J  xF (x)  dx 

5 

0.778 

0.423 

6 

0.704 

0.401 

7 

0.645 

0.380 

8 

0.596 

0.361 

9 

0.555 

0.344 

10 

0.520 

0.329 

11 

0.490 

0.315 

c 
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For  S  «  7  we  find  that 

xQ  =  0. 6  (3.30) 

Thus,  if  the  slender  body  load  is  truncated  at  xQ  =  0.6  (point  A  in  Figure 
3.1),  the  total  lift  is  the  same  as  that  obtained  with  the  present  theory  for 
Aq  =  1/2.  However,  the  equivalent  base  area  load  distribution  bears  very 
little  resemblance  to  measured  data  like  those  in  Figure  3.2. 

To  conclude  this  section,  we  calculate  a  typical  value  of  A^  from 
experimental  data.  The  lift  and  moment  coefficients  versus  angle  of  attack 
for  the  USS  Akron  Bare  Hull  and  complete  configuration  are  presented  in  Figure 
2.3  (also  from  Ref.  6).  The  coefficients  are  normalized  as  follows: 

c,  - 

q»(V)2/3 


where  RQ  is  the  maximum  body  radius.  We  compare  Bqs.  (2.11),  (3.24)  and 

(3.32)  to  obtain  the  following  expression  for  A^: 


dCN/da 


(3.33) 


From  Figure  3.3  we  have 


dC„ 

=  .0057/deg  -  .33/rad 


(3.34) 
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►  •.aaa.*iia 


ANGLE  OF  ATTACK 


Figure  3.3  Verietion  of  Lift  and  Pitching  Monent  Coefficients 
with  Angle  of  Attack  (Ref.  5) . 


and 


e  -  1/5.9  -  .169  (N  »  6) 


(3.35) 


Thus 


-  .45  (3.36) 

We  conclude  that  the  results  shown  in  Figure  3.1  for  A£  *  o.5  are  typical  for 
body  fineness  ratios  of  order  6  or  7. 

The  results  presented  in  this  report  have  focused  on  the  proper 
formulation  of  inviscid  slender  body  theory.  In  conclusion,  we  remark  that 
the  viscous  integral  equations  corresponding  to  our  fundamental  result, 
Eq.  (2.32),  can  be  written  down  at  once  following  the  procedure  outlined  in 
Section  II  E  of  Ref.  1.  We  simply  replace  the  Cauchy  singular  part  of  the 
kernel  function  with  its  viscous  form  to  get 


£n|x-y| 

. 


+  e 


o(x-y) 


K0(o|x-y| ) 


dy 


T  (x )  a  (x ) 

(3.37) 


where  KQ(z)  is  the  modified  Bessel  function  and  0  is  the  Reynods  number  based 
on  an  "appropriate"  body  length  scale  (e.g. ,  the  half  length  or  a  typical 
radius).  Hie  final  step  is  the  development  of  the  viscous  theory  in  order  to 
solve  Eq.  (3.37)  and  calculate  A^  as  a  function  of  Reynolds  number.  This  step 
will  be  a  breakthrough  insofar  as  our  fundamental  understanding  of  the  origin 
of  lift  on  a  body  of  revolution  and  in  our  ability  to  calculate  lift  from 
first  principles.  It  is  anticipated  that  the  viscous  work  can  be  completed  in 
a  future  research  effort. 


24 


IV.  CONCLUSIONS 


To  calculate  the  inviscid  lift  distribution  on  a  slender  body  of 
revolution,  the  essential  Cauchy  singularity  must  be  retained  in  the 
integral  equation  that  relates  the  load  and  upwash  functions. 

The  asymptotic  slender  body  theory  derived  herein  .has  an  exact  solution 
for  a  prolate  spheroidal  body.  The  eigenfunction  is  a  hypergeometric  or 
incomplete  Beta  function. 

The  lift  coefficient  for  a  slender  body  is  proportional  to  the  coefficient 
of  the  eigenfunction  in  complete  analogy  with  the  theory  of  lift  on  wing 
sections. 

The  coefficient  of  the  eigenfunction  is  of  order  unity  for  typical  slender 
bodies. 

The  qualitative  features  of  the  load  distribution  are  in  much  better 
agreement  with  experimental  data  than  the  distribution  obtained  with  the 
equivalent  base  area  concept  applied  to  conventional  slender  body  theory. 

The  theory  can  readily  be  modified  to  include  viscosity  following  the 
procedure  outlined  in  Section  II  E  of  Ref.  1. 
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